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Abstract - The passage of DNA through a nanopore can be effectively decomposed into two 
distinct phases, docking and actual translocation. In experiments each phase is characterized by 
a distinct current signature which allows the discrimination of the two events. However, at low 
voltages a clear distinction of the two phases is lost. By using numerical simulations we clarify 
how the current signature associated to the docking events depends on the applied voltage. The 
simulations show that at small voltage the DNA globule enhances the pore conductance due to an 
enrichment of charge carriers. At high voltage, the globule drains substantial charge carriers from 
the pore region, thereby reducing the overall conductance. The results provide a new interpretation 
to the experimental data on conductance and show how docking interferes with the translocation 
signal, of potential interest for sequencing applications. 


Introduction. — In the last decade, nanopore-based 
biosensing has become a burgeoning research field thanks 
to the impressive burst in the capability to fabricate de¬ 
vices based on solid state and biological pores [3||^. 
The working principle of the device is, in essence, simple. 
The nanopore connects two chambers containing an elec¬ 
trolyte solution. Under an electric voltage, ions migrate 
from one chamber to the other, with a conductance that 
depends on the details of the pore and the electrolytic 
solution. Since a translocating macromolecule alters the 
ionic flux, its passage can be detected and possibly the lo¬ 
cal nature of its monomers can be read off. Translocation 
of DNA is the most studied process being a promising 
technology for low cost/high throughput DNA analysis, its 
realization is still confronted with challenges [7 10. In 
parallel, a separate set of studies have focused on protein 
and polypeptide translocation, with the goal of structural 


characterization 11-15 


Over the pore extension, the DNA segment has a virtu¬ 
ally rod-like conformation and during the threading pro¬ 
cess the ionic current increases for low salt (< 0.4 M KCl) 
and decreases for high salt concentration [^[^ as com¬ 
pared to the open pore (so called free-pore in the following) 
case. The accepted explanation for the current reduction 
is that at high salinity the ionic motion is sterically ham¬ 
pered by the presence of DNA that reduces the effective 
pore section 18 . At low salinity instead, the counteri¬ 


ons cloud is more diffused and the charge carriers more 
mobile, an effect that prevails over the steric hinderance, 
increasing the conductance as compared to the free-pore 
condition 


16 17 


Recently, Kowalczyk and Dekker [^ and Vlassarev and 
Golovchenko 19 reported of a new type of hybrid event 


with two distinct signals, a current decrease followed by a 
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Figure 1: A) Sketch of the system. A cylindrical pore of length 
L and diameter d connects two reservoirs of diameter dr, with 
total length L^, = L + 2Lr, and aligned along the x axis. Bulk 
concentrations and the electric potential are imposed at the left 
(—) and right (+) boundaries. An immobile globule of total 
negative charge Q is located near the pore entrance. Note that 
the sketch does not reflect the actual aspect ratio of the system, 
see Fig. 


current increase. Interestingly, the amount of current re¬ 
duction is comparable to that of current enhancement and 
the transition between the two regimes takes place on a 
timescale shorter than the translocation event. The exper¬ 
iments were performed in low salt conditions and, hence, 
the increase in current was associated to the threading 
of DNA, while the decrease of current was interpreted as 
due to the presence of the DNA at the pore entrance. 
In fact, since the electric field is highly focused in the 
pore and access regions, the modulation of ionic current 
takes place only when DNA is in close proximity to the 
pore, in the so-called docked configuration. In addition, 
current reduction was attributed to the frictional forces 
exerted by the docked DNA on the solution [^. How¬ 
ever, the physical mechanism responsible for decreasing 
the ionic current remains elusive. In fact, while frictional 
forces can reduce electrokinetic transport, the transition 
between docked and threading DNA should be gradual. 
If friction has a leading role, a decrease of conductance 
should persist during DNA threading, since large portions 
of the DNA coil sit on both the cis and trans chambers of 
the device. 

Assessing the contribution to conductance stemming 
from a docked DNA is relevant to correctly discriminate 
the ionic signal arising from the translocating nucleotides. 
Due to the comparable intensity of ionic currents associ¬ 
ated to docked and translocating conformations, subtract¬ 
ing the systematic contribution originating from partially 
translocated DNA globules is particularly useful in the 
low voltage regime, when translocation is slow and there¬ 
fore the base pair content is amenable to detection. The 
present study aims at understanding the origin of cur¬ 
rent modulation by employing computer simulations in a 
model representation of the system. We consider an ideal¬ 
ized version of the device with DNA modelled as a charged 
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B 
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Figure 2: The current fields J", J"*' and J~, for the neutral, 
positive and negative species, on the y-z plane. A: J^, driven 
by electro-osmosis, on average directed ( -t), B: J+, on 

average directed { < -h) and C: JT, on average directed ( —>■ 

-t). Currents are reported in units of pbuik'^T 


globule docked near the pore entrance. 

Our results reproduce qualitatively the modulation of 
current at varying voltages as observed in experiments. 
In addition, they suggest an alternative cause for the ob¬ 
served modulation that is based entirely on electrokinetic 
effects. In particular, i) a substantial depletion of ions 
takes place within the pore due to the asymmetric charge 
drainage exerted by the globule, ii) the current decrease 
can be explained in terms of charge rearrangements, 
without the need for considering excluded volume or 
hydrodynamics interactions between DNA and the pore, 
iii) at low voltage the current increases as compared to 
the free-pore case, a regime that has not been detected 
by experiments. We further suggest that at low voltages 
the current increase observed in experiments should be 
assigned to hybrid docked/translocating configurations. 


System set-up and simulation method. — We 

consider a cylindrical nanopore of length L and diame¬ 
ter d — 0.3 L connecting two reservoirs of length = 3L 
and diameter d = L, (Fig. [^. The origin of the co¬ 
ordinate system is located at the center of the pore and 
the x-axis coincides with the pore axis. A potential dif¬ 
ference AV is applied between a positive and a negative 
electrode, placed at x = ±3.5 L, respectively. The fluid is 
composed of three species: positive charge carriers, neg¬ 
ative charge carriers and a neutral species 20,21 . The 


density fields are indicated as p~ and p” respectively. 
The electrolytic solution is assumed to be neutral and ho¬ 
mogeneous far away from the pore, with = p^uik 
the two open boundaries (x = ±3.5L). The Bjerrum and 
Debye lengths are set to Ib = 0.035 L and Xd = 0.05 L, 
respectively. These values are chosen in order to map the 
conditions employed in the experiments 

The role of the wall charge in the translocation exper¬ 
iment is two-fold: to enhance the electro-osmotic current 


^In Ref. L ~ 20 nm and d ~ Onm while. Is = 0.7 nm and 
lu = 10 nm, corresponding to a 0.15 M solution of a monovalent 
salt in water. 
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Figure 3: Positive (|/'''|, triangles) and negative (|/~|, inverted 
triangles) current intensities as functions of the applied voltage. 
Squares refer to the electrical current intensity (7c = |7^ —7“|). 
A non linear dependence of the electrical current intensity is 
apparent (7c ~ A17“ with a. = 0.93 ± 0.01). The current is 
represented in units of Ih = 7r(dr/2)^pj|"„;j,/uT. 



AV/mV 

Figure 4: Current in the presence of the globule at the pore 
entrance normalized by the value in absence of the globule, as 
a function of the applied voltage. Blue inverted triangle and 
red triangles refer to negative and positive carriers respectively, 
while green squares refer to the total electric current. 


by altering the amount of free charge carriers via the Don- 
nan effect and to introduce additional interactions with 
the macromolecule 15 22 . In the present work we ne¬ 


glect the presence of surface wall charge in order to keep 
the number of simulation parameters to a bare minimum 
and to isolate the effect of the docking globule on the ionic 
current. Such approximation allows capturing the leading 
effects in the modulation of the ionic conduction, while 
the presence of additional free charge carriers can be ac¬ 
counted for by modifying the bulk molarity. Secondly, as 
we are interested in docking conditions, the surface charge 
is expected to play a minor role on the DNA globule and 
charge redistribution outside the pore region. 

The presence of DNA docked at the pore entrance is 
mimicked by a spherical distribution of negative charges 
is positioned in a fixed location near the pore entrance, 
the total globule charge being Q (Fig. [^. 

The motion of the electrolytic solution is described by 
the density, velocity and electric fields of each species, hav¬ 
ing number densities = p",/9+,p“, equal unit mass, 
valence z® = 0,1,1, and equal mobilities. We adopt the 
continuum approach, implicitly assuming that the motion 
of individual molecules is not relevant for the problem 
at hand. Several studies indicate that for simple fluids 


this is a good description down to the nanoscale 23 24 


The electrokinetic flow is simulated by means of the Lat¬ 
tice Boltzmann Method, which reproduces the advection- 
diffusion dynamics for three species, a neutral one rep¬ 
resenting the solvent, and two charged species for the 
electrolytes, in the presence of electrostatic interactions 
20 21 25 . Simulations were performed with the parallel 


multiphysics software MUPHY 26 


Atomistic simulations and recent experiments indicate 
that the slip length for water on smooth surfaces hardly 
exceeds 1 nm 27 28 , therefore no-slip boundary condi¬ 


tions are applied for each fluid species at the solid wall 
27 28 . The materials usually employed for fabricating 


the nanopore are hydrophilic, thus allowing to safely disre¬ 
gard slip effects at the scale of interest. For electrostatics, 
the Poisson equation, V^(/> = — ^ (ep+ — ep~ + 0), where 
9 is the globule charge distribution, is solved with appro¬ 
priate boundary conditions. The tangential component of 
the electric field is continuous at the liquid-solid interface. 
Given the large disparity in the dielectric constants for 
water (e ~ 80) and the wall (e ~ 3.9 for silicon oxide and 
~ 7.9 for silicon nitride), we impose Neumann boundary 
conditions, V$ • n = 0. where n is the unit vector orthog¬ 
onal to the wall. Molecular Dynamics simulations have 
shown that the dielectric constant of water confined in a 
cavity of diameter 5nm is 71, close to the value in the 
bulk 29 . Hence, given the size of our system (d = 6nm), 


a minor alteration of the dielectric constant is expected 
and we employed e = 80 in the whole system. 


The docked globule is represented as a system of reg¬ 
ularized point-wise charges, each carrying a unit charge 
0.1 e, and distributed uniformly within a sphere of diame¬ 
ter d, i.e. the same diameter of the pore. The interaction 
strength is given by the Bjerrum length, Ib = e^/dTrefcsT, 
which is taken to be 0.7 nm, corresponding to liquid water 
conditions at ambient temperature. In order to resolve the 
double layer structure and avoid compressibility artifacts 
arising from the numerical method, we set the spacing 
of the computational grid to 0.02 nm. Fully developed 
flow is assumed at the inlet and at the outlet sections, 
each having a prescribed electrical potential (^i„ = 0 and 
4>out = AH). Dirichlet boundary conditions enforce the 
bulk densities of each species at the inlet and outlet, in 
particular, p+ = = p-^ = p-^^. 

For the Lattice Boltzmann method, the electrolytic 
solution is represented by three distribution functions. 
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f°‘{x,v,t), representing the probability of having a par¬ 
ticle of species a at position x, with velocity v and at time 
t. The species evolve according to the equation 


pz^ 3 

dtr +V-vr= a;(/“ - /“) + — • —/“ = C“[{/}] 

( 1 ) 

3/2 

/9“e is the Maxwell- 


where /“ = 


_ 1 _ 

2'7rv‘i 


Boltzmann local equilibrium, vt the thermal velocity, 
u = is the barycentric fluid velocity and the 

species velocity. The quantity ui controls the rate of re¬ 
laxation towards local equilibrium and is related to the 
kinematic viscosity i/ as v = v‘^{ljuj — 1/2). The distri¬ 
butions /“ are discretized in space over a Cartesian mesh 
and expanded in velocity space over a Hermite basis set. 
By Hermite projection, the distributions are replaced by 
f°‘{x,v,t) —>• /“ and analogously C°‘{x,v,t) —>• C“. The 
dynamics is generated according to the updating scheme 


fp{x + Cp,t+l)- f^{x,t) = C^{x,t) (2) 

where Cp represents a set of discrete speed that connects 
a mesh point to its neighbors. The species density is 
computed as p°‘ = J2p fp the species current as 

J“ = = Yhp^pfp- Example of the obtained cur¬ 

rent fields are shown in Fig. 

Results. — The current intensity of the species a is de¬ 
fined as /“ = p°‘u'^dA where A{x) is the local section 
perpendicular to the pore axis, p°‘ is the number density 
of a given species and the x-component of the velocity 
field. The total ionic current is Ic = —1~ ■ Asa prelim¬ 

inary check, we verified that in stationary flow conditions 
the system is globally neutral, and that in the absence 
of the globule, the number of positive and negative car¬ 
riers is exactly the same. Given the geometrical symme¬ 
try and the identical valence and mobility of the charged 
species the positive charge flow (electrodes (-I-) —)■ (—)) 
exactly balances the flow of negative carriers (electrodes 
(—) —>■ (-1-)), so that no electro-osmotic flow, i.e. the net 
motion of the neutral specie, is present. Moreover, the sys¬ 
tem response is found to be purely Ohmic and the ionic 
currents are individually proportional to the applied volt¬ 
age V (data not shown). In the following, we refer to 
quantities in the absence of the globule (free-pore case) 
with the subscript 0. 

The presence of the globule results in an abundance of 
positive carriers necessary to maintain the condition of 
global electroneutrality. The unbalance between positive 
and negative carriers reflects in a difference between \I'^\ 
and |/“|, see Figj^ 

To appreciate the role of the applied voltage, it is in¬ 
structive to compare the data in the presence and in the 
absence of the globule. Fig. [^reports the ratio between 
the current intensities for total globule charges Q = 100 e 
and Q = 0. In the explored voltage range the presence 
of the globule at the pore entrance leads to an increase 


of positive carrier current, > 1 (red triangles in 

Fig. 15 and to a depletion of the negative carrier one, 
I~ /Iq < 1 (blue inverted triangles in Fig. [ 5 . This mod¬ 
ulation is caused by the already mentioned abundance of 
positive carriers. 

A second observation is that both I~^/Iq and I~ /Iq 
decrease as the applied voltage AH increases, consistently 
with the depletion of carriers in the pore region as the 
voltage increases, see Fig. [5]4. Interestingly at low volt¬ 
ages Ic/Icfi > 0 , i.e. the presence of the charged globule 
results in an increase of the current, while at high voltage 
Ic/Ic,o < 0. The conductance modulation as a function 
of AH is interpreted by considering the simple case of a 
globule placed in the middle of a cylindrical tube of length 
Lx = L + 2Lr and diameter dr = L, i.e. a single cylinder 
with the same diameter of the reservoir with no interposed 
nanopore, spanning the whole length of the system. At 
zero voltage, the charge distribution is symmetric around 
the globule and, outside the Debye layer (i.e. at distance 
> 3 A_d from the globule), both and p~ recover the 
bulk value imposed at the two open boundaries. Under a 
finite voltage, the counterion cloud surrounding the glob¬ 
ule is asymmetrically distorted towards the negative elec¬ 
trode. The counterion density is therefore larger towards 
the negative electrode (charge surplus) than towards the 
positive electrode (charge depletion), see Fig. [^I. Clearly, 
outside the Debye layer the solution recovers local electro¬ 
neutrality {p'^ = p~), with an antisymmetric arrangement 
of charges with respect to the globule position induced by 
the applied voltage. 

The presence of the pore breaks the fore-aft symmetry 
and results in a substantial carrier depletion inside the 
pore (i.e. towards the positive electrode with respect to 
the globule). Such depletion is one order of magnitude 
more intense than the corresponding increase towards the 
negative electrode, see Fig. The pore induces an im¬ 
perfect electrostatic screening as compared on the oppo¬ 
site side of the globule. The effect increases with AH and 
hence, at high voltages, the carrier density inside the pore 
is strongly depleted by the charged globule. 

In essence, the docked globule produces two concomi¬ 
tant effects: i) the global increase in the available charge 
carriers and ii) the local depletion of carriers inside the 
pore. At low AH, the depleted pore plays a marginal role 
on the overall conduction and the leading effect is the in¬ 
crease of positive carriers, hence Ic > Ic,o- At high AH, 
the depleted pore dominates on the increase of the avail¬ 
able counterions, hence Ic < /c, fl¬ 
it is illustrative to relate the results of the present 
simulations with experimental findings of Kowalczyk and 
Dekker [^. We remind that these authors reported the 
intensity and the average duration of the current modula¬ 
tions, as shown in Fig. HP- The current decrease was asso¬ 
ciated to docking while the current increases to translo¬ 
cation events. Let us now focus on the current de¬ 
crease and the differential conductance, defined as AG = 
{Ic — /c.o)/AH. The experimental docking conductance 
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Figure 5: A: density profiles of positive (red triangles) and 
negative (blue inverted triangles) species normalized by the 
bulk value, for applied voltage AV = 5.2, 38.8, 77.5 mV (the 
arrows indicate the trend at increasing voltage). Inset: zoom- 
out of the profiles, with the light shaded area showing the data 
reported in the main graph and the dark shaded area the pore 
region. B: same as in the upper panel but in absence of the 
pore {AV = 18.1, 77.5mC). 



decreases monotonically in the range 200 mV < AV < 
000 mV (red squares in Fig. [6p) in qualitative agree¬ 
ment with our numerical result in the voltage range where 
Ic/Icfi <1- At small voltages {AV < 200mV), the ex¬ 
periments do not detect any current decrease, whereas the 
simulations show that Ic/Ic,o > 1- Moreover, experiments 
show an optimal enhancement at F ~ 125 mV (black cir¬ 
cles in Fig. 1^). The apparent discrepancy between ex¬ 
periments and simulations at low voltages can be recon¬ 
ciled according to the following argument. For 200 toF < 
AF < 600 mF the experiments detect hybrid events, with 
a clear-cut separation between reduction/docking (of du¬ 
ration td) and enhancement/translocation (of duration tt) 
events. Two distinct signals, called Id and It respectively, 
are measured as averages over the time intervals td and tt 
(Fig. [^). At low voltage, the simulations indicate that 
the ionic current is enhanced as compared to the free-pore 
case. Since both the docking and the translocation events 
feature an increase in conductance (Fig. |^) the docking 
events are hardly distinguishable from the translocation 
ones. In particular, as voltage diminishes, the docking 


events become slower than the translocating ones, as re¬ 
ported in [^. 

To estimate the conductance for low AF, we average the 
experimental data of over docking and translocating 
events, as 


AG = 


AGdtd + AGttt 
td + tt 


(3) 


where AGd, AGt, td and tt depend individually on the ap¬ 
plied voltage. We evaluate the conductance in the range 
50 toF < AF < 125 mF by using for td and U the expres¬ 
sions of ref. , and fit the experimental data in the inter¬ 
val 200 mF < AF < 600 mF to extrapolate the differen¬ 
tial conductances AGd and AGt at low voltage. In pass¬ 
ing, we note that the differential conductance due to dock¬ 
ing events is zero for AF ~ 130 mF. This value is com¬ 
parable to the numerical one (~ 30mF), notwithstanding 
the approximations involved, in particular regarding the 
large disparity of the globule size between simulation and 
experiments and the absence of pore wall charge. The fi¬ 
nal result of the averaging of eq. (1) is reported in Fig. 
1 ^ (blue solid line), which exhibits a good agreement with 
the experimental result, thereby lending good confidence 
in the proposed interpretation. 

In conclusion, the present work underscores the impor¬ 
tance of electrokinetic effects on DNA translocation exper¬ 
iments. In particular, it indicates that DNA docked near 
the pore entrance can either enhance or reduce the pore 
conductance at different applied voltages as compared to 
the free-pore case. Such non trivial phenomenon depends 
on the imperfect electrostatic screening that depletes the 
pore region from charge carriers. On the other hand, as 
the macromolecule begins threading, the conductance is 
dominated by the DNA thread in the pore region that in¬ 
creases and compensates for the charge drainage exerted 
by the DNA coil at the pore entrance. In the future, it will 
be interesting to investigate the role played by the charged 
wall and how the interplay between coil induced drainage 
and thread induced enrichment affects ionic conductance 
during DNA translocation. 


* * * 


Computing time at the CINECA supercomputing cen¬ 
ter (ISCRA grants FLEXPROT and NAPS) is kindly ac¬ 
knowledged. 


References 

[1] Russo C. J. and Golovchenko J., Proceedings of the 
National Academy of Sciences, 109 (2012) 5953. 

[2] Storm A., Chen J., Ling X., Zandbergen H. and 
Dekker C., Nature materials, 2 (2003) 537. 

[3] Hall A. R., Scott A., Rotem D., Mehta K. K., Bay- 
ley H. and Dekker C., Nature nanotechnology, 5 (2010) 
874. 

[4] Mohammad M., Iyer R., Howard K., McPike M., 
Borer P. and Movileanu L., Journal of the American 
Chemical Society, 134 (2012) 9521. 


P-5 






























M. Chinappi et al. 


A) low voltage 


Id [- ; 




1 - 1 It 

1 

lo td 

tt 

td 

lo 

!' 



B) high voltage 

It 


Id 


lo 



Voltage (mV) 


Figure 6: Upper panels: low (a) and high (b) voltage scenarios. 
At low voltage both docking and translocation events increase 
the current. Lower panel: experimental differential conduc¬ 
tance [^, associated to docking (red squares) and translocation 
(black circles). The dashed curves are fits of the experimental 
data for 200 mU < AV < 600 mU. The blue line is the new 
estimate, eq. for 50 mV < AV < 125 mV. 


Yang P., Nano Letters, 5 (2005) 1633. 

[18] Smeets R. M., Keyser U. F., Krapf D., Wu M.-Y., 
Dekker N. H. and Dekker C., Nano Letters, 6 (2006) 
89. 

[19] Vlassarev D. M. and Golovchenko J. A., Biophysical 
Journal, 103 (2012) 352. 

[20] Marini Bettolo Marconi U. and Melchionna S., 
Langmuir, 28 (2012) 13727. 

[21] Marconi U. M. B., Melchionna S. and Pagonabar- 
RAGA I., The Journal of Chemical Physics, 138 (2013) 
244107. 

[22] Pedone D., Firnkes M. and Rant U., Anal. Chem., 81 
(2009) 9689. 

[23] Chinappi M., Melghionna S., Casgiola C. and SuGCi 
S., The Journal of Chemical Physics, 129 (2008) 124717. 

[24] Bocquet L. and Charlaix E., Chemical Society Re¬ 
views, 39 (2010) 1073. 

[25] Melghionna S. and Marconi U. M. B., EPL (Euro¬ 
physics Letters), 95 (2011) 44002. 

[26] Bernaschi M., Melchionna S., Succi S., Fyta M., 
Kaxiras E. and Sircar J. K., Computer Physics Com¬ 
munications, 180 (2009) 1495. 

[27] Chinappi M. and Casciola C., Physics of Eluids, 22 
(2010) 042003. 

[28] Attard P., arXiv preprint arXiv:1302.6630, (2013) . 

[29] Senapati S. and Chandra A., The Journal of Physical 
Chemistry B, 105 (2001) 5106. 


[5] Laszlo a. H., Derrington I. M., Brinkerhoff 
H., Langford K. W., Nova I. C., Samson J. M., 
Bartlett J. J., Pavlenok M. and Gundlach J. H., 
Proceedings of the National Academy of Sciences, 110 
(2013) 18904. 

[6] Kasianowicz j. J., Brandin E., Branton D. and 
Deamer D. W., Proceedings of the National Academy 
of Sciences, 93 (1996) 13770. 

[7] Pennisi E., Science, 343 (2014) 829. 

[8] Kowalczyk S. W. and Dekker C., Nano letters, 12 
(2012) 4159. 

[9] Fyta M. G., Melchionna S., Kaxiras E. and Succi 
S., Multiscale Modeling & Simulation, 5 (2006) 1156. 

[10] Schneider G. and Dekker C., Nature biotechnology, 30 
(2012) 326. 

[11] Rodriguez-Larrea D. and Bayley H., Nature nan¬ 
otechnology, (2013) . 

[12] Bacci M., Chinappi M., Casciola C. and Cecconi E., 
The Journal of Physical Chemistry B, 116 (2012) 4255. 

[13] Bacci M., Chinappi M., Casciola C. M. and Cecconi 
F., Physical Review E, 88 (2013) 022712. 

[14] Madampage C., Tavassoly O., Christensen C., Ku- 
MARi M. and Lee J., Prion, 6 (2012) 110. 

[15] Cressiot B., Oukhaled A., Patriarche G., 
Pastoriza-Gallego M., Betton J., Auvray L., 
Muthukumar M., Bacri L. and Pelta J., ACS nano, 
6 (2012) 6236. 

[16] Chang H., Kosari F., Andreadakis G., Alam M., 
Vasmatzis G. and Bashir R., Nano Letters, 4 (2004) 
1551. 

[17] Fan R., Karnik R., Yue M., Li D., Majumdar A. and 


p-6 

























